Four Ways to Fit an Ion Channel Model

© 2019 Biophysical Society Mathematical models of ionic currents are used to study the electrophysiology of the heart, brain, gut, and several other organs. Increasingly, these models are being used predictively in the clinic, for example, to predict the risks and results of genetic mutations, pharmacological treatments, or surgical procedures. These safety-critical applications depend on accurate characterization of the underlying ionic currents. Four different methods can be found in the literature to fit voltage-sensitive ion channel models to whole-cell current measurements: method 1, fitting model equations directly to time-constant, steady-state, and I-V summary curves; method 2, fitting by comparing simulated versions of these summary curves to their experimental counterparts; method 3, fitting to the current traces themselves from a range of protocols; and method 4, fitting to a single current trace from a short and rapidly fluctuating voltage-clamp protocol. We compare these methods using a set of experiments in which hERG1a current was measured in nine Chinese hamster ovary cells. In each cell, the same sequence of fitting protocols was applied, as well as an independent validation protocol. We show that methods 3 and 4 provide the best predictions on the independent validation set and that short, rapidly fluctuating protocols like that used in method 4 can replace much longer conventional protocols without loss of predictive ability. Although data for method 2 are most readily available from the literature, we find it performs poorly compared to methods 3 and 4 both in accuracy of predictions and computational efficiency. Our results demonstrate how novel experimental and computational approaches can improve the quality of model predictions in safety-critical applications

Interpolatory methods for H∞\mathcal{H}_\infty model reduction of multi-input/multi-output systems

We develop here a computationally effective approach for producing high-quality $\mathcal{H}_\infty$-approximations to large scale linear dynamical systems having multiple inputs and multiple outputs (MIMO). We extend an approach for $\mathcal{H}_\infty$ model reduction introduced by Flagg, Beattie, and Gugercin for the single-input/single-output (SISO) setting, which combined ideas originating in interpolatory $\mathcal{H}_2$-optimal model reduction with complex Chebyshev approximation. Retaining this framework, our approach to the MIMO problem has its principal computational cost dominated by (sparse) linear solves, and so it can remain an effective strategy in many large-scale settings. We are able to avoid computationally demanding $\mathcal{H}_\infty$ norm calculations that are normally required to monitor progress within each optimization cycle through the use of "data-driven" rational approximations that are built upon previously computed function samples. Numerical examples are included that illustrate our approach. We produce high fidelity reduced models having consistently better $\mathcal{H}_\infty$ performance than models produced via balanced truncation; these models often are as good as (and occasionally better than) models produced using optimal Hankel norm approximation as well. In all cases considered, the method described here produces reduced models at far lower cost than is possible with either balanced truncation or optimal Hankel norm approximation

The relationships between bone mineral density in the spine, hip, distal femur and proximal tibia and medial minimum joint space width in the knees of healthy females

SummaryObjectiveTo investigate the relationships between bone mineral density (BMD) in the hip, spine, distal femur and proximal tibia and minimum joint space width (mJSW) in the knees of healthy women.MethodsWomen 22–68 years old without a history of knee pain, bone or joint disease or injury underwent a single, fixed-flexion knee X-ray. Radiographs were graded according to the Kellgren–Lawrence scale and analyzed for mJSW using a computer algorithm. Dual X-ray absorptiometry scans of the spine, hip, distal femur and proximal tibia were also acquired for each participant. Femur and tibia scans were acquired and analyzed using a modified version of the lumbar spine software.ResultsForty-five females, mean [standard deviation (SD)] age and body mass index (BMI) of 40.1 (13.9) years and 24.6 (4.5)kg/m2, respectively, participated. The mean (SD) mJSW was 4.64 (0.68)mm. Linear regression analyses controlling for age and BMI revealed that BMD in the femoral trochanter and the central two regions of the tibia (T2 and T3) was significantly related to mJSW in the knee. A backwards regression analysis performed to determine which region of interest is most significantly related to mJSW revealed that femoral trochanter BMD (β-value=0.416) is the most significant.ConclusionsIn contrast to the suggestion that BMD is negatively correlated with mJSW in the knees of osteoarthritic individuals, these results suggest that increasing BMD in the femoral trochanter and tibia is significantly associated with increasing mJSW in healthy females. Further investigation of this relationship is warranted

Considering discrepancy when calibrating a mechanistic electrophysiology model

Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions—that of uncertainty in the model structure or the equations themselves. The difference between imperfect models and reality is termed model discrepancy, and we are often uncertain as to the size and consequences of this discrepancy. Here, we provide two examples of the consequences of discrepancy when calibrating models at the ion channel and action potential scales. Furthermore, we attempt to account for this discrepancy when calibrating and validating an ion channel model using different methods, based on modelling the discrepancy using Gaussian processes and autoregressive-moving-average models, then highlight the advantages and shortcomings of each approach. Finally, suggestions and lines of enquiry for future work are provided. This article is part of the theme issue ‘Uncertainty quantification in cardiac and cardiovascular modelling and simulation’

Measurement of the View the tt production cross-section using eμ events with b-tagged jets in pp collisions at √s = 13 TeV with the ATLAS detector

This paper describes a measurement of the inclusive top quark pair production cross-section (σtt¯) with a data sample of 3.2 fb−1 of proton–proton collisions at a centre-of-mass energy of √s = 13 TeV, collected in 2015 by the ATLAS detector at the LHC. This measurement uses events with an opposite-charge electron–muon pair in the final state. Jets containing b-quarks are tagged using an algorithm based on track impact parameters and reconstructed secondary vertices. The numbers of events with exactly one and exactly two b-tagged jets are counted and used to determine simultaneously σtt¯ and the efficiency to reconstruct and b-tag a jet from a top quark decay, thereby minimising the associated systematic uncertainties. The cross-section is measured to be: σtt¯ = 818 ± 8 (stat) ± 27 (syst) ± 19 (lumi) ± 12 (beam) pb, where the four uncertainties arise from data statistics, experimental and theoretical systematic effects, the integrated luminosity and the LHC beam energy, giving a total relative uncertainty of 4.4%. The result is consistent with theoretical QCD calculations at next-to-next-to-leading order. A fiducial measurement corresponding to the experimental acceptance of the leptons is also presented

Search for strong gravity in multijet final states produced in pp collisions at √s=13 TeV using the ATLAS detector at the LHC

A search is conducted for new physics in multijet final states using 3.6 inverse femtobarns of data from proton-proton collisions at √s = 13TeV taken at the CERN Large Hadron Collider with the ATLAS detector. Events are selected containing at least three jets with scalar sum of jet transverse momenta (HT) greater than 1TeV. No excess is seen at large HT and limits are presented on new physics: models which produce final states containing at least three jets and having cross sections larger than 1.6 fb with HT > 5.8 TeV are excluded. Limits are also given in terms of new physics models of strong gravity that hypothesize additional space-time dimensions

Search for TeV-scale gravity signatures in high-mass final states with leptons and jets with the ATLAS detector at sqrt [ s ] = 13TeV

A search for physics beyond the Standard Model, in final states with at least one high transverse momentum charged lepton (electron or muon) and two additional high transverse momentum leptons or jets, is performed using 3.2 fb−1 of proton–proton collision data recorded by the ATLAS detector at the Large Hadron Collider in 2015 at √s = 13 TeV. The upper end of the distribution of the scalar sum of the transverse momenta of leptons and jets is sensitive to the production of high-mass objects. No excess of events beyond Standard Model predictions is observed. Exclusion limits are set for models of microscopic black holes with two to six extra dimensions

Measurement of the correlation between flow harmonics of different order in lead-lead collisions at √sNN = 2.76 TeV with the ATLAS detector

Correlations between the elliptic or triangular flow coefficients vm (m=2 or 3) and other flow harmonics vn (n=2 to 5) are measured using √sNN=2.76 TeV Pb+Pb collision data collected in 2010 by the ATLAS experiment at the LHC, corresponding to an integrated luminosity of 7 μb−1. The vm−vn correlations are measured in midrapidity as a function of centrality, and, for events within the same centrality interval, as a function of event ellipticity or triangularity defined in a forward rapidity region. For events within the same centrality interval, v3 is found to be anticorrelated with v2 and this anticorrelation is consistent with similar anticorrelations between the corresponding eccentricities, ε2 and ε3. However, it is observed that v4 increases strongly with v2, and v5 increases strongly with both v2 and v3. The trend and strength of the vm−vn correlations for n=4 and 5 are found to disagree with εm−εn correlations predicted by initial-geometry models. Instead, these correlations are found to be consistent with the combined effects of a linear contribution to vn and a nonlinear term that is a function of v22 or of v2v3, as predicted by hydrodynamic models. A simple two-component fit is used to separate these two contributions. The extracted linear and nonlinear contributions to v4 and v5 are found to be consistent with previously measured event-plane correlations

Search for dark matter produced in association with a hadronically decaying vector boson in pp collisions at sqrt (s) = 13 TeV with the ATLAS detector

A search is presented for dark matter produced in association with a hadronically decaying W or Z boson using 3.2 fb−1 of pp collisions at View the MathML sources=13 TeV recorded by the ATLAS detector at the Large Hadron Collider. Events with a hadronic jet compatible with a W or Z boson and with large missing transverse momentum are analysed. The data are consistent with the Standard Model predictions and are interpreted in terms of both an effective field theory and a simplified model containing dark matter

The performance of the jet trigger for the ATLAS detector during 2011 data taking

The performance of the jet trigger for the ATLAS detector at the LHC during the 2011 data taking period is described. During 2011 the LHC provided proton–proton collisions with a centre-of-mass energy of 7 TeV and heavy ion collisions with a 2.76 TeV per nucleon–nucleon collision energy. The ATLAS trigger is a three level system designed to reduce the rate of events from the 40 MHz nominal maximum bunch crossing rate to the approximate 400 Hz which can be written to offline storage. The ATLAS jet trigger is the primary means for the online selection of events containing jets. Events are accepted by the trigger if they contain one or more jets above some transverse energy threshold. During 2011 data taking the jet trigger was fully efficient for jets with transverse energy above 25 GeV for triggers seeded randomly at Level 1. For triggers which require a jet to be identified at each of the three trigger levels, full efficiency is reached for offline jets with transverse energy above 60 GeV. Jets reconstructed in the final trigger level and corresponding to offline jets with transverse energy greater than 60 GeV, are reconstructed with a resolution in transverse energy with respect to offline jets, of better than 4 % in the central region and better than 2.5 % in the forward direction